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Abstract 

We present a numerical simulation study of the density-dependence (p=2. 2-^4.0 g/cm 3 ) of the 
high-energy collective dynamics in vitreous silica at mesoscopic wavevectors (Q=l-^18 nm _1 ). The 
density-dependence of the longitudinal and transverse current spectra provides evidence that the 
excess modes observed in the density of states of this and many other glasses, i.e. the Boson Peak, 
arises from the high-Q limit of the quasi- transverse acoustic branch. This conclusion emerges from 
the comparison of the numerical results with the experimentally observed energy-shift and intensity 
variation of the Boson Peak with increasing density. 

PACS numbers: 61.43.Bn, 61.43.Fs, 63.50,+x 
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Insulating disordered solids, when compared to their crystalline counterparts, exhibit 
some common peculiarities in their low-temperature thermal properties and low-energy spec- 
troscopic features 1 ^, and in particular (i) a larger specific heat at temperatures up to ~1 
K, ascribed to tunneling processes^ 4 -, (ii) a much smaller thermal conductivity, which also 
shows a plateau in the temperature range ~l-=-10 K— , (iii) a quasi-elastic light- and neutron- 
scattering intensity, and, most important, (iv) an excess of modes in the vibrational density 
of states, known as the Boson Peak (BP). An unambiguous understanding of the origin 
of these extra modes, and of their possible relation with other reported anomalies, is still 
missing, in spite of the extensive research effort primed by the pioneering work of Buchenau 
et al.- and continued by many author o 6 ! 7 ^ 9 ! 10 ^ 1 ^ 12 ! 13 ^ 4 . No definite conclusion can yet be 
drawn as to the nature of the vibrational eigenvectors of the modes responsible for the 
BP: according to different authors, in fact, they are either spatially localized^, spatially 
delocalized and propagatin g 6 ' 9 ^ 3 , spatially delocalized but diffusive in character—. Even 
more important, still we do not completely understand why should disorder accumulate vi- 
brational eigenvalues in the same energy region, in such a broad variety of chemically and 
physically different materials. Some recent theoretical works have been devoted to this sub- 
ject, and among them we recall (i) the work of Elliott and coworkers, who assigned the BP 
in glasses to the lowest-energy van Hove singularity of the corresponding crysta l 14 i 16 ; (ii) the 
work of Grigera et al., who interpreted the BP as the precursor of the dynamical instability 
expected in a disordered structure as function of density^ 7 -; and (iii) the work of Gotze and 
Mayr and that of Schilling et al., who obtained a spectral feature recalling the BP within a 
mode-coupling-like description of the high frequency dynamics of a model glas o 1Sl19 . 

In this Letter, we present a Molecular Dynamics (MD) study of the high-frequency dy- 
namics in vitreous silica (v-Si02) at different densities in the range p=2.2-=-4.0 g/cm 3 . We 
demonstrate that the BP in vitreous silica arise from the high Q portion of the transverse 
acoustic branch. Our results quantitatively support and better specify i) the Elliott and 
coworkers' proposal 1 ^ 6 - that the BP is connected to the van Hove singularity of crystalline 
quartz; and ii) the Buchenau finding^ that the modes at the BP are local rotation of Si02 
tetrahedra. 

Experimentally, in Si02 at room pressure the BP appears as a broad peak in the plot 
of the Debye-normalized density of states g(E)/E 2 , centered at ^5 meV and having a 
width of ~ 5 meV (Fig. la). By increasing the sample density up to p=4.0 g/cm 3 , what is 
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observe d 20 i 21 i 22 i 23 is a high-energy shift of the peak and its concomitant intensity decrease 
(Fig. lb,c). In Fig. 1, the lines are experimental data, while symbols are the results of the 
present simulations (see below), and we would like to stress the excellent agreement between 
the two. It is important to realize that the excess of modes in the density of state itself 
g{E), corresponding to the peak at ~5 meV in g(E)/E 2 of Fig. 1, is actually centered (at 
room pressure) around « 15 meV, the shift being due to the 1/E 2 factor and to the very 
slow decay of the high frequency tail of the peak in Fig. 1 . With respect to this point, it 
is important to underline that -in the high Q limit and in the one-excitation (one-phonon) 
approximation- the dynamic structure factor S(Q,E) becomes proportional to g(E)/E 2 , 
and hence the longitudinal current spectrum Cl(Q, E)=E 2 /Q 2 S(Q, E) becomes oc g(E). 
Therefore, the modes giving rise to the BP must be searched, in the current spectra, around 
15 meV. This observation matches with the high Q panels of Fig. 2, where the computed 
currents (vide infra) show a broad bump at about this energy value. 

The investigated system consists of 680 SiC>2 units (iV=2040 ions), enclosed in cubic boxes 
of different lengths (from L=3.1359 nm, corresponding to p = 2.2 g/cm 3 for the glass at 
room pressure, down to L=2.5693 nm corresponding to a density of 4.0 g/cm 3 ), with periodic 
boundary conditions. The ions interact through the BKS 24 two-body interaction potential; 
the long-range interaction was treated by the Ewald sum technique. As it has already 
been demonstrated, this system reproduce quantitatively the high frequency dynamics of 
vitreous silica^. The glass configuration at room pressure was obtained by standard MD 
methods for lowering the temperature down to 300 K starting from a well equilibrated liquid 
configuration at T = 6000 K, followed by a conjugate gradient geometrical minimization on 
the potential-energy hypersurface for an accurate location of the minimum. Such minimum 
configuration was taken as the starting point to generate a series of compressed systems. 
At each compression step, the box size was scaled by ~1.5%, then the system was made to 
relax, and the new minimum configuration was searched by the conjugate gradient method. 
This procedure was repeated until the final density of 4.0 g/cm 3 , corresponding to a sample 
under an hydrostatic pressure of about 35 GPa 2 ^, was reached. A complete study of the 
structural and dynamical changes occurring during the compression will be reported in a 
forthcoming paper. We focus here on the changes in the high frequency dynamics that take 
place as the density is increased. 

The vibrational dynamics in the minimum configurations was computed in the harmonic 
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approximation by diagonalizing the dynamical matrix, to obtain the eigenvalues (E p ) and 
eigenvectors (e p (i)) of the p-th normal mode (p — 1 -j- 3N). From these quantities all 
vibrational characteristics can be derived. In particular, we have computed the density of 
states g(E) and the longitudinal and transverse currents spectra (C(Q,E)) which, in the 
one-excitation approximation, are given by: 

C2 (Q,E) = \K B Tl^MjAp\ Z P W£(Q)5(E - E p ) 

where 77 = L,T, a, (3 indicate Si and O, and W£(Q) is the spatial power spectrum of the 
(longitudinal or transverse) component of the eigenvectors: 

W p L (Q) = |S l (Q-e p «)exp(zQ-i? l )| 2 (1) 
Wp(Q) = |S J (gxe p (z))exp(2Q-i? l )| 2 . 

Here Q=Q/\Q\. 

The calculated g(E)/E 2 are reported in Fig. 1 together with the corresponding experimen- 
tal curves. The two sets of data compare favorably to each other, indicating the suitability 
of the employed potential model to follow the density dependence of the high frequency 
dynamics. In Fig. 2 we report longitudinal and transverse current spectra at selected Q 
values in the uncompressed sample (p=2.2 g/cm 3 ). For Q values larger than about 8 nm" 1 , 
both C L (Q, E) and C T (Q, E) show two distinct excitation maxima, a feature that becomes 
more and more evident at increasing Q. The excitation at higher energy disperses with Q 
and is observed at all Q values in the longitudinal current spectra, while it shows up as a 
weak shoulder in the transverse current spectra only at Q > 10 nm -1 . In agreement with 
previous finding o^ 2 ^ 2 ? , we assign this feature to the longitudinal sound-like branch. The 
behaviour of the low-energy excitation is in some sense complementary: it is always present 
in the transverse current spectra, while it appears in the longitudinal current spectra only 
at Q > 8 nm" 1 . At small Q, the low-energy peak disperses with a sound velocity of ~ 3800 
m/s (appropriate for the transverse sound mode), and becomes almost non- dispersing at 
Q > 8 nm -1 (Fig. 3). We will call this low-energy feature -which is the main feature in the 
transverse current spectra- the transverse acoustic mode. 

The presence of the signature of transverse dynamics in the longitudinal current spectra, 
and vice-versa, is only apparently surprising. Indeed, the polarization character of the modes 
(which is better and better defined at increasing wavelength, i. e. when the vibration sees 
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the medium as an elastic continuum) becomes ill denned at short wavelengths (mixing phe- 
nomenon). This is at the basis of the growth of peaks associated to the opposite polarization 
modes in the current spectra, and of the increased visibility of these peaks at increasing Q 
values. In Fig. 3 we report the computed dispersion curves of transverse excitations (open 
symbols) for samples of three densities (p = 2.2, 2.7, 4.0 g/cm 3 ), and of longitudinal excita- 
tions (crossed circles) of the uncompressed system^. In the same figure (full symbols), are 
also reported -at the available Q values, the experimental excitation energies (maxima of 
the current spectra) of the uncompressed Si0 2 glass, as measured by IXS and INS 28 i 30 . The 
agreement between MD and experimental data clearly indicates that the peaks observed 
in the INS experiment at Q larger than 8 nm" 1 - , must be associated with the transverse 
dynamics, which appears in the measured (longitudinal) spectra due to the mixing. 

In Fig. 3, we also observe a strong positive dispersion of the velocity of the (longitudinal) 
sound waves, which continue to energy well above the BP energy. This disper- 

sion, observable also in MD simulations of Lennard- Jones glasses^ and in previous v-SiC>2 
calculations^, can be due to the interaction of the sound waves with a relaxation process^ 
or with other modes. The investigation of the origin of the positive dispersion of vl, and is 
possible interplay with the BP modes, is beyond the purpose of the present work. 

For our present purposes, the first important result emerging from Fig. 3 is that the 
"transverse" branches (open symbols) in low- and intermediate-density samples, at large Q's, 
flatten at high Q's to an energy value which increases with increasing density (~15 meV at 
room condition and Q=15 nm 4 ). The density of vibrational states (g(E)), associated with 
the branches which flatten, will have an excess of modes with respect to the Debye behaviour 
at these energies, reminiscent of the van Hove singularity of the corresponding crystal. As 
mentioned, the corresponding peak in g(E)/E 2 , will be red-shifted at lower energies (~5 
meV at room condition). Therefore, we conclude that the Boson peak originates from the 
modes associated to the flat portion of the acoustic transverse dispersion curve. 

This result is in agreement with the recent theoretical work of S. Taraskin et al., who 
associate the BP to the glassy counterpart of the lowest-energy van Hove singularity of the 
corresponding crystalline structure^. In this respect, it is worth noting that the transverse 
acoustic branch at Q larger than ^8 nm -1 is the glassy counterpart of a transverse optic 
phonon branch of a-quartz (almost flat at ^4 THz, i. e. ~16 meV). This branch, in the 
extended Brillouin zone scheme which is more appropriate for disordered materials, is the 
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prosecution of the transverse acoustic branch 3 - 2 -. According to Boysen at al.— , this branch is 
(close to the M-point) strongly temperature-dependent, and its softening is responsible for 
the a-/3 transition in quartz. More importantly, the atomic displacements induced by the 
lattice modes of this branch in quartz, as determined in 3 - 2 -, are very similar to the frustrated 
localized rotation of SiC>4 tetrahedra, that in vitreous silica -according to Buchenau et al£- 
are the modes contributing to the BP. 

The assignment of the BP to the flattening of the transverse acoustic branch is strength- 
ened by the behaviour of the high-energy transverse dynamics upon densification. Spectra 
similar to those of Fig. 2, but for the sample at the highest studied density (p=4.0 g/cm 3 ), 
are reported in Fig. 4, while the corresponding dispersion curves are reported in Fig. 3. One 
can observe that (i) the two current spectra are now "pure", and no evidence of wrong- 
polarized modes is present; (ii) the positive dispersion of the L-branch is absent; and, more 
importantly, (iii) the T-branch no longer shows a flattening. This latter observation, to- 
gether with the experimental and numerical finding that the BP intensity strongly decreases 
on increasing densit y 21122123 (as evident in Fig. 1), gives a decisive support to the finding 
that the BP is produced by the flattening of the quasi- TA branch. 

In conclusion, by comparing the spectra of the longitudinal and transverse current, and 
by studying their density dependence, we have shown that the BP in v-SiC>2 is to be ascribed 
to the quasi-transverse acoustic modes, whose dispersion relation becomes Q-independent 
at high Q, and consequently, gives rise to an excess of modes with respect to the Debye 
behaviour. The whole picture presented here reconciles different previous studies on the 
origin of the BP: i) the high Q part of the transverse acoustic branch in vitreous silica 
-which gives rise to the BP- is the counterpart of a low lying transverse optic branch in 
a-quarta 3 - 2 .; ii) according to Boysen at al.— the softening of this branch at the M-point 
give rise to the a-to-(3 transition; iii) in agreement with i) and ii), according to Taraskin 
et dlM^&, the BP arises from the softening of the lowest-energy van Hove singularity of 
the corresponding crystals; iv) according to Dorner and coworkers the eigenvectors of this 
branch in quartz correspond to rotation of SiC>4 tetrahedra; v) finally, in agreement with 
i) and iv), according to Buchenau et al.— the modes contributing to the BP are localized 
rotation of Si0 4 tetrahedra. 

This work was supported by INFM Iniziativa di Calcolo Parallelo, and by MURST Pro- 
getto di Ricerca di Interesse Nazionale. One of us (GR) greatly acknowledge illuminating 
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discussions with U. Buchenau and B. Dorner. 
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FIGURE CAPTIONS 



FIG. 1: Normalized density of states g(E)/E 2 of vitreous silica at three different densities 
p: (a) 2.2 g/cm 3 (full line and circles); (b) 2.7 g/cm 3 (dashed line and squares); (c) 4.0 g/cm 3 
(triangles). Lines are experimental data from Ref.— , symbols are the results of the present 
simulation. 

FIG. 2: Selected examples of longitudinal (full lines) and transverse (dashed lines) current 
spectra at the indicated Q-values in nm _1 , for the sample at p=2.2 g/cm 3 . 

FIG. 3: Dispersion relation of the main peaks appearing in the current spectra. Crossed 
circles are the maxima of the longitudinal current at p=2.2 g/cm 3 (the line is a guide for 
the eyes). The open symbols refer to the maxima of the transverse current at p=4.0 g/cm 3 
(down-triangles), p=2.7 g/cm 3 (squares), p=2.2 g/cm 3 (up-triangles). Full symbols are INS 
(squares) and IXS (circles and diamonds) experimental data for p=2.2 g/cm 3 taken from 
Ref.s^ao. 

FIG. 4: Selected examples of longitudinal (full lines) and transverse (dashed lines) current 
spectra at the indicated Q- values (in nm" 1 ) for the highest density sample (p=4.0 g/cm 3 ). 
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Figure 2. Pilla et al. 
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Figure 3. Pilla et al. 
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Figure 4. Pilla et al. 
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